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ABSTRACT 

The origin of intergalactic magnetic fields is still a mystery and several scenar- 
ios have been proposed so far: among them, primordial phase transitions, structure 
formation shocks and galactic outflows. In this work we investigate how efficiently 
galactic winds can provide an intense and widespread "seed" magnetisation. This may 
be used to explain the magnetic fields observed today in clusters of galaxies and in the 
intergalactic medium (IGM). We use semi-analytic simulations of magnetised galactic 
winds coupled to high resolution N-body simulations of structure formation to esti- 
mate lower and upper limits for the fraction of the IGM which can be magnetised up 
to a specified level. We find that galactic winds are able to seed a substantial fraction 
of the cosmic volume with magnetic fields. Most regions affected by winds have mag- 
netic fields in the range 10~ 12 < B < 1CU 8 G, while higher seed fields can be obtained 
only rarely and in close proximity to wind-blowing galaxies. These seed fields are suf- 
ficiently intense for a moderately efficient turbulent dynamo to amplify them to the 
observed values. The volume filling factor of the magnetised regions strongly depends 
on the efficiency of winds to load mass from the ambient medium. However, winds 
never completely fill the whole Universe and pristine gas can be found in cosmic voids 
and regions unaffected by feedback even at z — 0. This means that, in principle, there 
might be the possibility to probe the existence of primordial magnetic fields in such 
regions. 
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1 INTRODUCTION 

The investigation of magnetic fields in the intergalactic 
medium (IGM) and in the intracluster medium (ICM) is 
a new challenge for modern cosmology. While magnetic 
fields have been successfully detected in clusters of galax- 
ies using a variety of techniques, such as Faraday rotation 
measurements and observations of diffuse large-scale syn- 
chrotron emission ( f or recent reviews see | Carilli fc Taylor! 
120021: IWidrowlEool iGovoni fe Ferettill2004) . only a hand- 
ful of pioneering observations suggest that magnetic fields 
of significant intensity may be present in the IGM 1 and no 
accurate measurements are available for those regions of the 
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1 We define the intracluster medium (ICM) as the gas bound to 

clusters and groups of galaxies, while we call intergalactic medium 

(IGM) the diffuse gas in the Universe which is not bound to any 

such structure. This distinction is clear in numerical simulations, 

but it may not always be so in observations. 



Universe with densities of order of the mean cosmic density 
or lower. 

Cluster magnetic fields detected through Farada y rota- 
tion measurements are of the order of a few uG. iKim et all 
yjMj) measured field strengths of about 2 u.G on scales 
of 10 kpc in a statistical sample of point sources ob- 
served through the Coma cluster, and Clarke, Kronberg 
and Bohringer (2001) derived field strengths of 4 to 8 
11G in a sam ple of point sources observed in 16 low red- 
shift clusters. Ijohnston-Hollitt fc Ekera J2004) have done 
a similar analysis of southern galaxy clusters and pro- 
vided an in dependent con firmation of the statistical signal 
reported bv IClarke et alJ <l2001f) . Magnetic field strengths 
of 10 to 40 11G have been reported for the cool cores 
of t he Hydra A and Cygnus A clusters on scales of 3- 5 
kpc (iDreher. Carilli fc Perlevlll987l iTavlor fc Perlevlll993l) . 
IVoet fc En filin (2005) apply a maximum likelihood estima- 
tor for magnetic power spectra, based on the theory of tur- 
bulen t Faraday screens jEnfflin fc VoetfeOOStlVogt fc Enfilir] 
2003), to a high quality Faraday rotation map of the 
north lobe of the radio source Hydra A produced by the 



Bertone et al. 



novel PACERMAN algorithm foolae. Voet fc Enfilir] 120051: 
IVoet. Dolae fe Enfilinl [20051) . They find a cluster central 
magnetic field strength of 7 ± 2uG. 

There have been few atte mpts to det e ct int ergalactic 
magnetic fields beyond clusters. Kim et alJ 111989]) detected 
radio-synchrotron radiation from an extended magnetic field 
in the region of the Coma supercluster, with a field strength 
of 0.1 to 0.01 uG. This radio emission comes from a re- 
gion with an enhanced number density of galaxies, which 
might indicate a group of galaxies in the process of merg- 
ing with Coma. iBagchi et al] 1120021) report possible radio- 
synchrotron evidence of intergalactic magnetic fields in the 
region ZwCL 2341.1+0000 at z ~ 0.3, where galaxies seem 
to be aligned along a filament. However, in both cases it re- 
mains controversial whether these measurements prove the 
magnetisation of the diffuse IGM or instead the state of the 
gas belonging to the supercluster environment. 

Theoretically, three basic mechanisms could explain the 
pres ence of magnet ic fields in the ICM and IGM (for a re- 
view |Widrowj|2002|) : i) they could have a primordial origin, 
ii) they were seeded in the IGM by battery effects of shock 
waves or inhomogeneous radiation fields, or iii) their exis- 
tence is linked to the formation and evolution of galaxies. A 
primordial origin implies that weak seed fields of about ~ 
10~ 18 G were prod uced during the early phas es of the cos- 
mic evolution fsee lReeall987l : IKro nbcrg 1994). These seed 
fields must have later been amplified by extremely efficient 
mechanisms, if the fields we are to observe today were of 
the order of a fraction of a fi.G. Alternatively, seed magnetic 
fields may have been produced in protogalaxies and galaxies 
and subsequently ejected into the ICM and IGM. In both 
cases, cosmic shear and/or dynamo effects may intervene to 
amplify the seed fields. 

Cosmological magnetohydrodynamic simulations of 
magnetic fields in gala xy clus t ers h ave been performed 
by iDolag. Bartelmann fc Leschl (11993) , assuming a homo- 
geneous seed field of 10 -9 G at z ~ 15. They find that 
seed fields of this intensity can be efficiently amplified to re- 
produce the magnetic field strength observed in local clus- 
ters. However, inhomogeneous seed fields can similarly lead 
to the right amount of cluster magnetisation: iMiniati et al] 
(2001a) follow the evolution of the magnetic fields produced 
by Biermann-battery effects in shock waves from the cos- 
mic structure formation process and find that they can be 
further amplified by cosmic shear flows, leading to strong 
cluster fields. 

Two different mechanisms can explain the transport of 
the magnetic seed fields produced in galaxies into the ICM 
and the diffuse IGM. The magnetic energy may be eas- 
ily ejected into the IGM and ICM by je ts and radio lobes 
emerging from powerful radio galaxies if ricjvjg] ^ffij: [R^gsJ 



1987; Dalv & Locb 1990; Chakrab arti. Rosner fc Vainshtein 
1994). Enfilin et al. (1997) estimate d the injection of mag- 
netic fields by radio galaxies, and IVolk fc Atovanl (2000) 
followed the subsequent amplification through dynamo ef- 
fects, finding that it can account quite accurately for the 
observed magnetic field strengths and cosmic ray densities in 
the ICM. Kronberg et al. (2001) make a step further and es- 
timate the magnetic flux transported directly into the IGM. 
The second mechanism is related to the ejection of magne- 
tised interstellar gas and supern ova ejecta by galactic winds. 
iKronberg. Lesch fc Hood (1199911 consider such a scenario of 



magnetised outflows from star forming galaxies and find that 
a substantial fraction of the intergalactic medium could have 
been seeded at relatively high redshift. However, their esti- 
mate was based on a power-law extrapolation of the cosmic 
star formation rate to high redshifts, which is not consistent 
with the known peak of star format ion activity at z ~ 1 
JMadau et alll99d) . |Pe Yound dl992l) follows the amplifica- 
tion of a variety of magnetic seed fields by turbulent dynamo, 
driven by the motion of galaxies in the ICM, and concludes 
that this mechanism cannot account for the magnetic fields 
observed in clusters. 

Galactic winds emerging from star forming galax- 
ies have been detect ed thro u ghout t he history of 
the Unive rse (e.g. iPhillips! Il993l, JH eckman et al] 
2000L ICecil. Bland-Hawthorn fc Veflleuxl l200l 

Frve Broadhurst fc Benitezj I2002J lAdelberger et al] 



20031 IVeilleux. Cecil fc Bland-Hawthornl |2005j) and are 

often advocate d to explain the metal enrichme nt of 

the IGM (e.g. lAguirre et all 1200 it iTheuns et all 120021 
iBertone. Stoehr fc White! 12003. BSW05 hereafter). Al- 
though it is still controversial if winds from large galaxies 
can fully account for t he observed level of m etal enrichment 
of the IGM at 2 ~ 3 dAguirre et al]|2005t) . it is clear that 
they do play a role in galaxy formation and there is no 
doubt whether they eject a consistent fraction of the metal 
and gas content of galaxies at any redshift. 

Magnetic fields are clearly observed in the interstellar 
medium of galaxies (Beck 2001), where they are believed to 
play an essential role during star formation. In a few cases, 
magnetic fields are also observed in the galact ic material 
outflowing from star forming galaxies ( Reuter e^flJ_[W94|: 
iGolla fc Hummell994lTullmann et all200CD . lReuter et^l] 
( 1992) estimate a magnetic field st rength of 10 /xG in th e ma- 
terial outflowing from M82, while IPahlem et al. (1997J) find 
a field strength of about 7 /xG in the outflow of NGC 4666. 
Birk, Wicchcn & Otto ( 1999) argue that galactic winds are 
able not only to advect galactic magnetic fields but also to 
amplify them by shear flows from Kelvin-Helmholtz insta- 
bility. If a significant fraction of the outflows are magnetised, 
it is reasonable to expect that they may at least partially 
contribute to the magnetisation of the ICM and IGM. The 
transport of the magnetic flux by galactic outflows may oc- 
cur at virtually any redshift: at higher redshifts the main 
contribution may come from protogalaxies with total masses 
of the order of > 10 6 Mq , while at lower redshifts galaxies 
with stellar masses comparable to the ones of local dwarf 
galaxies or larger may become the dominant source of mag- 
netic fields. 

In this paper, we use high resolution simulations of 
galaxy formation, which include semi-analytic prescriptions 
for the physics of galactic winds, to investigate the efficiency 
of outflows to seed the ICM and IGM with magnetic fields 
and cosmic rays. Numerical models of galactic winds do par- 
tially explain the transport of the metals produced in galax- 
ies into the IGM and, since there is an analogy between the 
transport of metals and the transport of magnetic fields, 
we investigate the possibility that galactic winds can eject 
seed magnetic fields from galaxies and transport them into 
the IGM and ICM. However, we do not intend to investi- 
gate in detail all the physical phenomena connected to the 
transport of the magnetic fields in our simulations, because 
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semi-analytic models are intrinsically not suitable to treat 
the magneto-hydrodynamics involved in full detail. 

We present two different scenarios for the ejection and 
transport of magnetic fields by galactic winds: a "conserva- 
tive" scenario and an "optimistic" scenario. For each case, 
we then analyse the effect of the efficiency of winds to af- 
fect the results. In the first scenario, we aim to provide a 
"conservative" estimate of the global magnetisation of the 
Universe and we only consider the ejection and the trans- 
port of magnetic "seed" fields by winds, while amplification 
mechanisms are neglected. As a consequence, these results 
should be treated as lower limits. In the second scenario, we 
include a simple prescription for the amplification of the seed 
magnetic fields during the wind evolution and we predict a 
level for the magnetisation of the Universe that represents 
an extreme case and that can be therefore regarded as an 
upper limit. 

Given the resolution of our simulations, we focus on the 
evolution of winds emerging from galaxies with M* > 10 8 
Mq at z < 10. We are interested in determining the fraction 
in volume of the ICM and IGM which can be magnetised by 
winds and we calculate the minimum magnetic energy and 
the minimum cosmic ray energy density which can be ejected 
into the ICM and IGM by winds. 

The paper is organised as follows. In Section [5] we de- 
scribe our set of simulations and how we implement our 
recipes for the transport and amplification of the seed mag- 
netic fields on top of the prescriptions for the evolution of 
galactic winds. In Section [3] we present our results for the 
transport of magnetic fields by winds and we discuss the 
possible amplification of the seed fields by cosmic shear and 
turbulence. We briefly discuss the ejection, transport and 
amplification of cosmic rays in Section|3]and finally we draw 
our conclusions in Section [H] 



2 THE SIMULATIONS 

In this work, we use the simulations of galactic winds pre- 
sented in BSW05 and we implement new semi-analytic 
recipes for the transport of the magnetic fields on top of the 
pre-existing scheme. We assume a ACDM cosmology with 
matter density Q, m — 0.3, dark energy density £7a = 0.7, 
Hubble constant h = 0.7, primordial spectral index n — 1 
and normalisation as — 0.9. 

In Subsection 12.11 we briefly describe the set of high 
resolution N-body simulations and the semi-analytic model 
for galaxy formation. We present our prescriptions for the 
physics of galactic winds in Subsection 12.21 and for the evo- 
lution of magnetic fields and cosmic rays in Subsections 12.31 
and 14. II respectively. 



2.1 Formation and Evolution of Galaxies 



A set of high resolution N-body simulations, "M3" llStoehrl 
2003), is used to follow the evolution of the dark matter 
and the formation of structures with time. The "field" sim- 
ulation M3 describes an approximately spherical region of 
space with diameter 52 /i _1 Mpc and average density close 
to the cosmic mean. The fact that only about half the en- 
closed galaxies at z — are field galaxies, while the rest are 
in groups and poor clusters, makes it particularly suitable 



to investigate the effects of feedback on the diffuse IGM. 
The M3 simulation efficiently combines a high resolution in 
mass with a large simulated volume, necessary to study the 
effects of feedback from galaxies with stellar masses as low 
as dwarfs in their proper cosmological context. The num- 
ber of particles in the simulated region is about 7 ■ 10 7 and 

the particle mass is 1.7 • 10 8 /i _1 Mm. The simulations were 

i 1 1 i w 

performed by Stochr ( 2003) using the parallel treccodc GAD- 
GET I iSpringel. Yoshida fc Whiteil200ll) . The dark matter 
evolution was followed from redshift z — 120 down to red- 
shift z — and 52 simulation outputs were stored between 
z = 20 and z = 0. 

The formation and evolution of galaxies is mod- 
elled with the semi- analytic technique proposed by 
Kauffman n et al.l 1119 99') in the new implementation by 
Springel et al.l (2001). Merging trees extracted from the sim- 
ulations are used to follow the galaxy population in time, 
while simple prescriptions for gas cooling, star formation and 
galaxy merging model the processes involving the baryonic 
component of the galaxies. At z = 3 a total of about four 
hundred thousand galaxies are identified and about three 
hundred and fifty thousand are present at Z — 0. The two 
largest clusters have a total mass of about 10 14 /i _1 Mq. 



2.2 Prescriptions for Galactic Winds 

BSW05 implement new prescriptions for the long-term evo- 
lution of winds in the semi-analytic code of ISpringel et alJ 
(2001). Their model aims to provide simple predictions 
about the effects of winds on the diffuse IGM and, in partic- 
ular, on its metal enrichment history. Here, we give a brief 
qualitative description of the model, but we refer the inter- 
ested reader to BSW05 for a more detailed discussion. 

The dynamics of galactic winds is modelled as a two- 
phase process: a pressure-driven, adiabatic phase, followed 
by a momentum-driven phase. We treat these two phases 
separately, according to the different physical processes that 
drive the expansion of winds. During the first phase, an 
over-pressurised bubble of hot gas emerges from a star form- 
ing galaxy and evolves adiabatically until its cooling time 
becomes shorter than its dynamical expansion time. This 
Sedov-Taylor phase of the evolution terminates when the 
loss of energy by radiation becomes substantial and the ma- 
terial at the outer edge of the bubble starts to cool down. 
During the second phase of the wind evolution, a thin shell 
of cooled material accumulates at the shock front and the 
dynamics of the wind becomes dominated by the momentum 
imparted by the outflowing material onto the thin shell. This 
second phase is the classic momentum-driven snowplough 
and both the shell and the outflowing wind material are 
cool. For simplicity, we make the assumption of spherical 
symmetry. 

According to lOstriker fc McKea l|l988t) . during the 
Sedov-Taylor phase of adiabatic expansion the dynamics of 
the outflow is described by the equation for the conservation 
of energy of a spherical bubble with energy injection at the 
origin: 

dE w 1 • 2 . a rfl 
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where_E w is the wind energy, R and v B are the radius and 
the velocity of the shock, M w and i> w the mass outflow rate 
and the outflow velocity of the wind, p , Po, v and u the 
density, the pressure, the outward velocity and the inter- 
nal energy of the surrounding medium, respectively. Mu is 
the total mass internal to the shock radius. The entrain- 
ment fraction e is a parameter that defines the fraction of 
mass that the bubble sweeps up while crossing the ambient 
medium. This equation states that the bubbles are powered 
by the material outflowing from the galaxies and that they 
are slowed down by gravity and by the ram pressure and the 
thermal pressure of the ambient medium. 

The equation of motion for a spherically symmet- 
ric thin shell is given b y the conservation of momentum 
( Ostriker fc McKeelllQSsi : 



d 
d-t {mV * 
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drt 

eA-kR p v (v s — V ), 



(2) 



with m, R and v B the mass, the radius and the velocity of 
the shell and (j> the gravitational potential of the dark matter 
halo. The first term on the right-hand side represents the 
momentum injected by the starburst, the second term takes 
into account the gravitational attraction of the dark matter 
halo and the two final terms represent the thermal and the 
ram pressure of the surrounding medium. 

The conservation of mass law for a pressure driven bub- 
ble is 
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while for a thin shell 
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Initial conditions for the mass loss rate and the velocity 
of winds are taken from Sh u. Mo Ik Mad (2005) and cal- 
ibrated to match the observations of winds both in local 
starbursts and in high redshift ga laxies (e.g. lMartinlll999l 
iFrve. Broadhurst fc Benitea I2002J) . M w and u w are calcu- 
lated as a function of the star formation rate of the galaxy 
blowing the wind: 



Ah 



M° 71 K M 



© y r 



= 320 • Ml 



'K 



-1/2 



km s 



(5) 



(6) 



where K is a constant that takes into account various prop- 
erties of the interstellar medium (ISM) like, e.g., the effi- 
ciency of conduction relative to the thermal conductivity 
of clouds, the minimum radius of clouds in the ISM and 
the dimension of star-forming regions (see IShu. Mo fc Mad 
2005 for a comprehensive discussion) . With this assumption, 
the momentum input M v v w is only weakly dependent on K 
(oc K 1 ' 2 ), with a stronger dependence on the star formation 



rate (oc M" 



. The energy input is proportional to the star 



formation rate, but does not depend on K and therefore on 
all other galaxy properties. 

Both K and e define the mass loading of winds, but 
while K determines the amount of gas entrained from the 
ISM of the galaxy, e is the fraction of gas swept up by winds 
from the ambient medium after blow out. The remaining 



fraction of mass 1 — e is assumed to be in dense clouds which 
are unperturbed by the wind. A low value e < 1 may reflect 
either a clumpy ambient medium or a heavily fragmented 
wind, while e ~ 1 describes a near-homogeneous medium, 
which can be entirely swept up by the wind. For this work 
we assume K — 0.5 and we consider two values of the en- 
trainment fraction parameter, that is e = 0.1 and e = 0.3. 
The first value represents the behaviour of winds with ineffi- 
cient mass loading, which can easily escape from galaxies at 
any redshift, while the second represents highly mass loaded 
winds, whose efficiency to expand far into the IGM is sig- 
nificantly reduced with respect to less mass loaded winds. 
BSW05 find that the results for winds only weakly depend 
on K, while the strongest dependence in the model is on the 
entrainment fraction parameter e. Currently, there are no 
observational constraints on the parameter e and our simu- 
lations are unable to con strain it on the basis of our results. 
iBertone fc White! (J2005) try to constrain the parameter by 
comparing the statistical properties of synthetic Lya spec- 
tra with observed ones, but are unable to put any constraint 
on either of the two models: both scenarios are equally likely 
and none is in contrast with observations. The volume fill- 
ing factor of galactic winds is similarly not accurately con- 
strained by observations: published estimates vary between 
0.004% and 40% and in principle neither of the models we 
present in this work can be excluded on the basis of obser- 
vations. 

When winds expand inside the haloes of galaxy groups, 
they sweep up the hot and possibly metal enriched intra- 
group medium and they are subject to the cluster gravita- 
tional attraction, thermal and ram pressure forces. Winds 
expanding into the IGM entrain metal-free gas, which is 
assumed to have constant density and pressure and whose 
velocity field is determined by the Hubble flow. In princi- 
ple, winds have to be far more energetic to escape from 
large dark matter haloes than to escape from field galax- 
ies or small groups, since the forces they have to overcome 
are stronger. In general, winds have a higher probability to 
escape the haloes of isolated galaxies than to escape from 
groups, and it is unlikely that they can escape from large 
clusters, although we cannot prove this in our simulations 
because of the lack of massive dark matter haloes. Ram pres- 
sure plays an important role to quench outflows and its in- 
fluence is particularly strong in large haloes, where it can 
prevent winds from blowing (Kaofcrcr et al. 2006). Thermal 
pressure in low-redshift haloes can also slow down winds 
considerably, but is normally negligible with respect to the 
other forces at high redshifts. 

2.3 Prescriptions for Magnetic Fields 

The aim of our calculations is to describe the injection of 
seed magnetic fields into the IGM by galactic winds. With 
our current means, this can only be done in an approxi- 
mative way. One of the main simplifications we use in the 
following is the assumption that the magnetic fields in winds 
are dynamically negligible. This does not imply that organ- 
ised magnetic fields at the base of the wind could not play 
a role in launching the outflows. However, we assume that 
the conversion of magnetic forces into kinetic energy hap- 
pens before the fields enter the expanding wind bubble. The 
remaining magnetisation evolves passively within the out- 
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flow, allowing a simplified description of its evolution, as 
we briefly describe below. A more detailed derivation of our 
equations is given in Appendix lAl 

We define the injection of the magnetic fields ejected by 
galaxies into the IGM as the magnetic field seeding, since this 
expression literally describes the process well. Potentially, 
shear flows in the IGM can subsequently amplify these seed 
fields to the observed levels. We note that the term mag- 
netic field seeding is often used in the literature to describe 
the Battery-effects, which generate magnetic fields without 
requiring a pre-existing seed magnetization, due to correc- 
tions to the ideal magneto-hydrodynamic equations. This 
initial generation of magnetic fields is believed to take place 
in stars, in galaxies and during the early phases of the evo- 
lution of the universe. In this work, we do not attempt to 
model the process, but we assume that it has happened in 
the past, allowing galaxies to create their own magnetic field 
and to amplify it to the observed level. This existing mag- 
netic field is then ejected by the galaxy and advected in to 
the IGM by a galactic-scale outflow. 

In our conservative scenario, we treat the seed magnetic 
fields analogously to the metals and we do not process them 
any further: this implies that we describe the injection of the 
seed fields by galaxies into the outflows and their dilution 
due to the expansion of winds, but that we do not take into 
account any possible mechanism that could produce ampli- 
fication. With this assumption, we can make predictions for 
the injection and transport of seed magnetic fields into the 
IGM and ICM, which effectively result in a lower limit of 
the actual magnetisation of the Universe. 

We present our prescriptions for the amplification of 
the seed fields in Subsection 12.3.21 where we describe the 
build-up of our optimistic scenario. We discuss other pos- 
sible mechanisms that could further amplify the seed fields 
transported by winds in Subsection 13.31 and we make a few 
consideration about the transport of cosmic rays in winds 
in Section 2] 

2.3.1 Conservative Scenario 

Since in this work we want to make order of magnitude pre- 
dictions, we neglect the possible influence of the magnetic 
fields on the evolution of the galactic winds and we set up 
the new semi-analytic prescriptions for the magnetic fields 
on top of the pre-existing scheme of BSW05. However, this 
approximation turns out to be well justified: in fact, in our 
calculations we find that the magnetic energy is always sev- 
eral orders of magnitude smaller than the total energy in 
the winds. 

To simulate the transport of magnetic fields by galactic 
winds, we need two main ingredients: an equation for the 
evolution of the magnetic fields and suitable initial condi- 
tions. To be consistent with the evolution of winds set by 
equations Q and J3J, we derived an equation for the evolu- 
tion of the magnetic energy in winds: 
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(7) 



Eb is the total magnetic energy in a bubble, -Es in is the rate 
of magnetic energy injected by winds, V w is the wind volume 
and Vw is the variation in the wind volume due to the ex- 
pansion. The second term on the right hand side of equation 



10 expresses the decrease in the magnetic energy due to the 
adiabatic expansion of the bubbles and is a negative contri- 
bution. This can be deduced from the induction equation for 
a frozen-in magnetic field expanding isotropically, when the 
shear amplification term is neglected. A detailed derivation, 
which discusses the underlying assumptions, can be found 
in Appendix^ Here we give a simplified argument. The in- 
duction equation yields B oc 1/R 2 and Eb oc V B 2 oc 1/R. 
This implies that 



(8) 



where the numerical factor 1/3 reflects the assumption of no 
special alignment between the local field direction and the 
local direction of the expansion of the wind. 

We assume that the magnetic energy injection rate is 
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where Vi n is the volume of the injected material, pi n is the 
density of the gas injected into the wind, £s in is the mag- 
netic energy density of the wind material outflowing from 
the galaxy, and the mass outflow rate M w is given by equa- 
tion ©. Equation|5|assumes that the magnetic field injected 
at the base of the wind has the strength of a typical galactic 
magnetic field. Since the production of a galactic-scale mag- 
netic field by, for example, an a — Q dynamo mechanism is 
not instantaneous, we assume that in our simulations mag- 
netic energy can be injected in winds only after the host 
galaxy has completed a full rotation. The rotational period 
of the galaxy is calculated as a function of the disk radius 
and circular velocity, that is i rot = 2tvR s /v c . 

In our model, winds make a transition between the pres- 
sure driven phase and the momentum driven phase when 
their cooling time becomes shorter than their dynamical 
timescale. However, when simulating the evolution of the 
magnetic fields in the wind, we continue to use equation 
Q for the conservation of the magnetic flux throughout the 
wind evolution. This is a reasonable approximation, because 
the magnetic energy can be injected into the winds indepen- 
dently of the dynamical phase of the wind and equation J7J 
does not critically depend on any quantity linked to a par- 
ticular phase of the wind. 

A crucial quantity in Eq. @ is the density pi n of the 
gas outflowing from a galaxy at a rate M w . A estimate of 
this quantity can be derived as a function of the mass loss 
rate M w of a wind through the entire surface area iirRg of 
a galaxy at the velocity « w : 
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where we assume that the galactic radius is a fixed fraction 
of the virial radius, that is R g — R200/W, as m BSW05. 
The magnetic energy density es in associated with the mass 
injected in winds can be expressed as 

4/3 
n, , - — I ^ I , (11) 



with Bgai the galactic magnetic field and pism the density 
of the interstellar medium (ISM) of the galaxy. Here we as- 
sume that the magnetic field lines expand adiabatically and 
isotropically when the ISM is blown out of the disk of a 
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galaxy and is injected into a wind. The factor in paranthcses 
in Eq. dl IB represents therefore the dilution of the magnetic 
fields when the density of the ejected gas drops below the 
density of the ISM. The adiabtic index 4/3 derives from the 
relationship B oc 1/R 2 in the induction equation, which im- 
plies B oc p 2 ' 3 and B 2 oc p 4 ' 3 for isotropic expansion. The 
density of the ISM is given by 

3Mism no s 

PISM - l^jp (12) 

with Mism the mass of the ISM. 

Galaxies in t he local Universe have magnetic fields o f 
order 5 ^G (e.g. iDahlem et alJll997t IChvzv et al"ll2003t) . 
while interacting galaxies can have magnetic fields as high 
as 30 uG (e.g. ISoida et alJ 1200 ll IChvzv fc Beckl I2004L 
lHarnett et alJ l2004l. In our simulations, we will express 



-Bgai = b ■ B, 



(13) 



where B — 1 p.G is a conservative reference magnetic field 
value at the base of the winds and 6 is a "fudge" factor to 
rescale our results to less conservative choices of an initial 
magnetic field. 

2.3.2 Optimistic scenario 

In this Subsection we describe our simple and purpose- 
fully optimistic prescription for the amplification of mag- 
netic fields during the wind evolution. To set up this pre- 
scription, we consider only the amplification caused by two 
mechanisms directly linked to the wind evolution, namely 
the injection of the shocked interstellar gas into the wind 
and the expansion of the wind into the surroundings of the 
source galaxy. 

The first mechanism we consider is the amplification 
of the magnetic fields at the base of the winds, where the 
outflowing material escapes the visible regions of galaxies. 
Here, the magnetic fields embedded in the ejected gas can be 
amplified through Kelvin-Helmholtz instabilities. This com- 
plex p rocess is described in detail bv lBirk. Wiechen fc Ottol 
( 1999) and we do not intend to disc uss it further. What is 
releva nt for this work is the fact that lBirk. Wiechen fc Ottol 
(1999), when simulating the process, find that the magnetic 
energy density increases by a factor of 3, almost indepen- 
dently of their set of initial conditions. We take this into 
account in our simulations by scaling the magnetic field we 
inject in winds by the fudge factor b — y/3 in Eq. I|13|l . In 
addition, we do know that observations of magnetic fie lds in 
local galaxies yield strengths of about -Bgai ~ 6|jG fe.g. lBeckl 
l200ll: IDahlem et alJll997t iGhvzv et aljfeooaft . For our opti- 
mistic scenario, we therefore set the magnetic field strength 
associated to the galactic ISM to be B ga i = 6uG. If we incor- 
porate this value of a more realistic galactic magnetic field 
in the fudge factor, together with the amplification factor of 
\/3 at the base of the wind by Kelvin-Helmholtz instability, 
we obtain a total fudge factor b — 10, which we include in 
our simulations. 

The second mechanism we consider for the amplifica- 
tion of magnetic fields is the effect of shear flows inside wind 
cavities and bubbles during the wind expansion into the am- 
bient medium. To estimate how strongly the magnetic fields 
could be amplified, we first need to determine the character- 
istic timescale of this mechanism. Since motions inside winds 



can take place on scales of the order of the wind radius R, 
an effective timescale r c fr can be estimated as 



-1 Us 

Teff = Ti- 



ll 



=ns-*w 



(14) 



where v B is the proper wind expansion velocity, vn = H (t)-R 
is the Hubble flow at distance R from the source galaxy and 
H (t) is the Hubble constant as a function of time. The term 
v s — vh represents the peculiar velocity of a wind with respect 
to the Hubble flow and is the only shear motion whose con- 
tribution to the amplification of the magnetic field in wind 
we consider, r) is a field amplification efficiency factor which 
depends on the poorly known flow geometry inside the wind. 
This amplification by shear motions directly counteracts the 
dilution of the magnetic fields due to the adiabatic expan- 
sion term and can be incorporated in Eq.|7|for the evolution 
of the magnetic energy as a positive term: 



-77 Eb = Eb u 
at 



Tott O Vw 



(15) 



If we assume that t] — 1 and we remember that the ratio 
Vw/ (3VW) can be expressed as v s /R, Eq.^] simplifies to 
d 



<lt 



Eb = Es in — H it) Eb- 



(16) 



In this extremely optimistic case in which r\ = 1 , the di- 
lution of the magnetic fields by adiabatic expansion and the 
amplification by shear motions cancel out, leaving behind 
only the term H (t) Eb- This factor is small with respect 
to the magnetic energy injection term -Es in and it therefore 
represents only a small correction in Eq. 1161 



3 RESULTS 

In this Section we show our results for the magnetisation of 
the IGM and ICM. In Subsection 13. II we discuss the simu- 
lated properties of the magnetic fields in winds for the con- 
servative and the optimistic scenarios. In Subsection 13. 21 we 
present our results for the fraction of magnetised IGM as 
a function of cosmic time and we qualitatively describe the 
effect of galaxy clustering on the distribution of the mag- 
netic fields in our simulated region. We present our results 
in terms of comoving magnetic field strengths, as it is usu- 
ally done in the community. Since our simulations produce 
physical estimates of the magnetic fields, we use the relation 
-Bcomoving = -Bphysicai / ( 1 + z) 2 to convert from one system of 
reference to the other. 

In the following Subsections, all results are presented 
for a total of four different models, which we list here for 
convenience. The first two models represent our conservative 
scenario and correspond to two different wind efficiencies: 
the "s = 0.1" model represents the case of efficient winds 
and the "e = 0.3" model the case of less efficient winds. 
There are then two additional models which represent our 
optimistic scenario. In the following, they will be referred 
to as the "e = 0.1 A" model and the "e = 0.3A" model and 
they correspond to efficient and less efficient wind models 
respectively. 

3.1 Magnetic Fields and Cosmic Rays in Winds 

Before describing our main results for the magnetisation of 
the IGM, we focus briefly on the properties of the magnetic 
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Figure 1. The magnetic field strength B w for the conservative scenario in winds as a function of the wind density for three different 
redshifts (z = 6, z = 1 and z = 0, from left to right) and for our two wind models with e = 0.1 (upper panels) and e = 0.3 (lower panels) 
in the conservative scenario. 



fields in winds as they emerge from our simulations. Our 
semi-analytic model computes the total magnetic energy in 
a wind Eb , which can easily be transformed in the root mean 
squared magnetic field strength _B W by the equation: 



Eb 



B 



4-7T 3 



(17) 



with R the wind radius. B„ should be regarded as a root 
mean square magnetic field strength within the wind rather 
than as the value of a uniform magnetic field. 

We find that B w depends on several factors: among 
them, the wind density and the total magnetic energy in- 
jected by the starburst. Consistently, we never observe val- 
ues of B larger than 1 u.G in our conservative scenario, which 
means that the dilution of the field lines due to the expan- 
sion of the wind is a significant effect at any stage of the 
wind evolution. 

In Fig. Q we use the two wind models in our conserva- 
tive scenario to show the comoving magnetic field strength 
in individual winds as a function of the wind density. Re- 
sults are presented for three different redshifts (z = 6, z = 1 
and z = 0, from left to right) and for our two wind models 
with e = 0.1 (upper panels) and e = 0.3 (lower panels). Each 
point represents a wind at the particular snapshot time. The 
total number of winds is between about 10 4 and about 3T0 5 , 
depending on the redshift and on the model. There is a clear 



trend that indicates that B w increases with the wind den- 
sity. The scatter is mostly due to the different stages of evo- 
lution of each wind: most winds are currently expanding into 
the IGM, but a number of them are recollapsing onto the 
galaxies, a phenomenon known as galactic fountain. Winds 
outflowing from galaxies which have only recently assembled 
most of their mass may not be magnetised, since the galaxy 
may not have had time to build up a magnetic field (see 
Subsection l2.3H . In Fig.0 the scattered points on the right- 
hand side of the distributions of magnetic fields indicate 
newly formed bubbles. These winds, which are still expand- 
ing in proximity of galaxies, normally have higher densities 
and their magnetic field strengths _B W depend more sensi- 
tively on the magnetisation of the material outflowing from 
the ISM of the galaxy. Young winds may have low B w when 
little magnetic energy is injected by the starburst because 
the wind mass loss rate is small. However, in a few cases, 
when the wind mass loss rate is substantial, they may also be 
highly magnetised (see the points in the upper right corner 
of the distribution). The points corresponding to the lowest 
wind densities belong either to weakly magnetised winds or, 
in most cases, to winds which have joined the Hubble flow. 
By changing the fudge factor b in the expression for B ga i, 
the results in Fig. can be rescaled by the same factor 6. 

In Fig. [5] we show the volume averaged RMS magnetic 
field strength in the Universe B n (left y axis) and the 



8 Bertone et al. 



o 



en 
O 



-1 



/ 


.. 


8 


^■^' .••'' - 






















9 


,' .••'' ^^^^ - " ~ : 




• jX***"^ "* 


















n 


- •'' /^ • *" p = D 1 - 




/" ' ' 




i X ,' e = 0.3 




s 


1 


: / £ = 0.1 A -. 




e = 0.3 A ; 


? 


, , : 



t„ 



(Gyr) 



-16-T 



18 ?> 

CD 



-20 £ 
s= 

CO 
en 

-22 o 



-24 



Figure 2. The volume averaged RMS comoving magnetic field 
(left axis) and the cosmic ray energy density (right axis, see Sub- 
section 14.11 in winds as a function of time for the two different 
scenarios and the two values of our wind model parameter e. If 
the input parameter B In is modified by a fudge factor b, the re- 
sults presented here for the cosmic RMS magnetic field strength 
can be rescaled by the same factor b, while the magnetic energy 
density can be rescaled by a factor b 2 . 



volume averaged cosmic ray energy density ecR (right y axis) 
as a function of the cosmic time. The cosmic ray energy 
density has the same behaviour as magnetic field strength 
because we assume the magnetic and the cosmic ray energy 
densities to be in equipartition. Results are presented for 
the optimistic and the conservative scenarios and for our 
two wind models. 

The volume averaged RMS value of the cosmic mag- 
netic fields increases with time in both scenarios and for 
both wind models, in agreement with the trend we observe 
in Fig.^for the mean value of the magnetic fields in individ- 
ual winds. As expected, the magnetic field strength increases 
more strongly in the optimistic scenario than in the conser- 
vative scenario: this is a consequence of the fact that in this 
scenario the dilution of the magnetic fields in winds because 
of the adiabatic expansion is largely compensated by the 
internal shear amplification. 

Fig. [5] also shows that the level of magnetisation of 
the Universe depends critically on the efficiency of winds to 
transport matter, cosmic rays and magnetic energy out and 
far from galaxies. As stated in Subsection 12.21 the winds in 
models with e = 0.1 are much more efficient in blowing out 
of galaxies than in other models. As a consequence, the level 
of magnetisation of the Univese produced in this model is 
more than an order of magnitude larger than in the e = 0.3 
models. 



3.2 Volume Filling Factors of Magnetic Fields 

In this Subsection we present our results for the volume fill- 
ing factors of magnetised regions. By volume filling factor we 
mean the fraction of our simulated region affected by winds 
with a desired value of the magnetic field strength. It is im- 
portant to estimate the filling factor of magnetised regions 
in this way, because it gives us the possibility to make sim- 
ple predictions about where the magnetic fields may have 
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Figure 3. The volume filling factors of magnetised regions as 
a function of the cosmic time. The four panels show the filling 
factors of regions with magnetic field strengths in different ranges 
for all four models. Results also indicate the filling factors of the 
cosmic ray energy density in winds (see Subsection 14. II . 



been ejected and transported by galactic winds and what 
kind of intensities we may expect them to have. Eventually, 
this gives us information about the ultimate mechanism that 
generated the magnetic fields we currently observe in the 
Universe. In the following, we will try to answer the ques- 
tions: can galactic winds eject the seed magnetic fields that 
have generated the magnetic fields in clusters and in the 
IGM? Or do we need other mechanisms to explain them? 

The volume filling factor /b of magnetised regions is cal- 
culated using the same technique used by BSW05 to obtain 
the volume filling factor of galactic winds. Indeed, for each of 
our two models, the sum of the filling factors of regions with 
different magnetisation intensities shown in Fig.|S|is exactly 
the volume filling factor f v of winds for the two models with 
(K;e) = (0.5; 0.1) and (K;e) = (0.5; 0.3) shown in Fig. 8 
of BSW05. To calculate /b we superpose a grid with 512 
pixels over our simulated region and we identify all those 
pixels which correspond to the desired magnetic field inten- 
sity. The volume filling factor is then simply given by the 
ratio between the number of pixels with the required mag- 
netic field strength B and the total number of pixels in the 
grid. The magnetic field strength in each pixel is recovered 
by summing up the magnetic energy density of each wind in- 
tersecting that particular pixel and afterwards converting it 
into B using equation ITU . This conversion assumes that the 
magnetic field energy is uniformly distributed inside winds: 
this is only statistically correct. 

Our estimates for the volume filling factors of magne- 
tised regions as a function of the cosmic time are shown in 
Fig. for all models. The four panels represent the volume 
filling factors of regions with magnetic fields in four differ- 
ent intervals: i) B > 10" 7 G; ii) 10" 8 < B < 10" 7 G, hi) 
10" 9 < B < 10" 8 G and iv) 10" 12 < B < 10~ 9 G, respec- 
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Figure 4. The relative volume filling factors of magnetised re- 
gions as a function of the cosmic time. The four panels show the 
filling factors of regions with magnetic field strengths in different 
ranges for all four models. Results also indicate the relative filling 
factors of the cosmic ray energy density in winds (see Subsection 

Era . 



tively. Magnetic fields with strengths higher than 10~ G 
have extremely low or low filling factors at all epochs in both 
scenarios. On the other hand, the total filling factor of mag- 
netic fields with strengths in the range 10 -1 < B < 10~ 
G is already significant at z ~ 3 and approaches unity at 
the present epoch in models with e = 0.1. This implies that 
in the e = 0.1 and in the e = 0.1 A models corresponding to 
efficient winds, the ejected magnetic fields permeate most of 
the universe at z ~ 0. 

In Fig. [I] we show the relative volume filling factors of 
magnetic fields with given strengths as a function of cosmic 
time. We calculate the relative volume filling factor of mag- 
netised regions with strengths in the same intervals used in 
Fig. |!J] This quantity represents the ratio between the filling 
factors /b shown in Fig.[3]and the total volume filling factor 
of galactic winds f v as a function of time, that is / = fb/fv 
/ v is the fraction of our simulated region affected by winds, 
independently of their magnetisation. Physically, / indicates 
the fraction of space affected by winds in which there is a 
magnetic field of a given intensity. It appears that in the 
conservative scenario only a tiny fraction (less than 10~ 4 ) of 
the volume filled by winds contains magnetic fields higher 
than 10 -7 G, while most of the wind filled regions have 
10" 12 < B < 10~ 8 G throughout the history of the Uni- 
verse. In the conservative models, most magnetised regions 
have a magnetic fields smaller than about 10 -9 G and only 
at the present time most regions can reach level of magneti- 
sations as high as B ~ 10~ 8 G. 

On the other hand, in the optimistic models, most winds 
have magnetic field strengths in the range 10 -9 < B < 10 -7 . 
In this models, at z — most magnetised regions appear to 
have a magnetic field in the range 10 -8 < B < 10 -7 G 



and almost no regions are only weakly magnetised (that is 
B < 10~ 9 G). In particular, in the e = 0.1 and in the e — 
0.1A models, most of the Universe is highly magnetised at 
z ~ 0. If this were the case, galactic winds would be proven 
to be an efficient mechanism to expl ain the magnetis ation 
of filaments and the IGM observed bv lKim et al.l lll989l) and 
iBaechi et al. (2002), if cosmic shear or other amplification 
mechanisms could further amplify the seed fields by a factor 
of 10-1000. 

Fig. |S|shows slices of the magnetic field evolution in our 
simulated region as a function of redshift. Results are pre- 
sented for the e — 0.3 (left panels) and e = 0.1 (right pan- 
els) models of the conservative scenario. In the optimistic 
scenario models with e = 0.3j4 and e — 0.1 A, results are 
qualitatively similar, but the colour table would have to be 
stretched to higher values of the magnetic field intensity. 
From top to bottom and from left to right, the redshifts of 
the snapshots are z — 3, z — 1, z — 0.5 and z = 0. We cut a 
thin slice with a side of 52 ft -1 Mpc through our simulated 
region and we colour-code the magnetic field strength from 
white (no field) to black (strongest fields). The value of the 
magnetic field in each pixel is the integrated value of the 
magnetic field B recovered in the 3D grid with 512 3 cells 
used to calculate the filling factors of magnetised regions in 
Subsection l3.2l There is a substantial difference between the 
maps of the magnetic fields in the two wind models shown 
here. Although the total volume filling factors / v differ by 
less than an order of magnitude, the effect is clearly visible 
in Fig. This again is a consequence of the efficiency of 
galaxies to blow winds, as determined by the parameter e. 
In the e — 0.1 model, winds fill almost all the simulated 
region, leaving us only a small chance, if any, to observe pri- 
mordial magnetic fields in the Universe at z ~ 0. On the 
contrary, models with higher entrainment fraction parame- 
ters and less efficient winds as the e = 0.3 model, show that 
although the winds do fill a large portion of the simulated 
volume, at least half or more of the intergalactic space is 
not affected by feedback from galactic outflows. This means 
that, in principle and with the appropriate experiment, one 
could observe primordial magnetic fields even in the present- 
day Universe. 

The distribution of the magnetic fields in Fig. fs] shows 
clearly that they are not statistically homogeneously dis- 
tributed in space, but that the magnetisation is somewhat 
correlated with the underlying large scale structure of mat- 
ter and galaxies. This is easily understandable, because most 
of the galaxies blowing winds form in regions of relatively 
high density, like filaments and groups, while fewer objects 
populate very low density regions like voids. Therefore, it 
may be possible that the plasma in voids is still pristine, 
allowing to search for uncontaminated primordial magnetic 
fields. The probability to detect pristine gas in low density 
regions at z = ranges from about 10% for the two models 
with e = 0.1 to about 70% for the two models with e = 0.3. 



3.3 Further amplification of the Seed Magnetic 
Fields 

In our simulations, most magnetised winds emerge from 
galaxies in filaments or galaxies that reside in proximity 
of regions with enhanced densities. Incidentally, these re- 
gions are also those where one would expect the amplifica- 
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Figure 5. The evolution of the magnetic field strength with cosmic time. Results are presented for the wind models with s = 0.3 (left 
panels) and with e = 0.1 (right panels) of our conservative scenario. The initial injected magnetic field is B lu = 1 u.G. In each figure, 
from top to bottom and left to right, the redshifts of the shown snapshots are z = 3, z = 1, z = 0.5 and 2 = respectively. The colour 
coding is the same for all snapshots and for both models. The comoving dimension of the simulated region is 52 h Mpc. 



tion mechanisms such as cosmic shear and dynamo effects 
to be more efficient, because of the stirring of shear flows 
and turbulence by the infall of matter on the newly forming 
structures. 

Clusters of galaxies are expected to be only weakly 
rotating, if at all, and other field amplification mech- 
anisms than the a — Q dynamo have to be found 
to explain the observed magnetic fields . iJaffd |l980) 
and iRuzmaikin. Sokolov fc Shukurovl 1(1989) suggested that 
magnetic fields could be amplified by the small scale turbu- 
lence prod uced by th e motio n of g alaxies through the ICM. 
Howev er, iDe Youna J1992T) and iGoldshmidt fc Rephaelil 
( 1993) demonstrate that this mechanism alone cannot fully 
produce the levels of magnetic fields observed in clusters. 

There is some observatio nal evidence fo r the exis- 
tence of turbulence in clusters. ISchuecker et alj (120041) find 
hints of the presence of turbulence in the Coma cluster by 
analysing the pressure fluctuati ons revealed by X-ray ob- 
servations. lEnBlin fc Vogtl (|2005l) suggest that hydrodynam- 
ical turbulence may be responsible for the generation of the 
Kolmogoroff-like magnetic power spec trum observed in the 
Hydra A cluster JVoet fc Enfflinll200FJ) . 

The merging of galaxy clusters is prob ably the most 
important source of turbulence i n the ICM iTri bblc 1993, 
Roettieer. Burns fc Stone I1999L iNorman fc Brvanl Il999). 
Roettieer. Stone fc Burns lll999l) find that magnetic fields 



can be amplified th ree times more e fficiently in cl usters with 
ongoing merge rs. Dolag, Bartel mann fc Leschl (1999) and 
iDolag. Bartelmann fc Leschl J2002l) simulate the magnetic 
field evolution in clusters and find that the field strength 
can be increased by factors of order 10 3 by compression and 
shear of random flows during the cluster formation. Such 
an amplification factor would be sufficient to bring even the 
seed magnetisation of a nG we find in our conservative mod- 
els to the observed (J.G level of magnetisation in galaxy clus- 
ters. 

Galactic winds emerging from galaxies in clusters can 



contribute to the magnetisation of the ICM in clusters even 
when they do not escape the gravitational attraction of the 
halo. In fact, when star formation occurs, some material can 
be shocked and ejected outside the plane of the galaxy. This 
gas is then quickly stripped away from the galaxy by the ram 
pressure of the hot cluster gas JSoker, Bregman fc Sarazinl 
Il99ll iGaetz. Salpeter fc ShavivMl987l) and mixed into the 
ICM. The stripped material would therefore contribute to 
chemically enrich, shock heat and magnetise the ICM of the 
cluster, although it would be unable to power a wind and 
reach the IGM. Unfortunately, we are currently unable to 
reproduce this effect in our numerical model. 

Outside clusters and bound structures, magnetic fields 
in winds can be amplified by turbulence, if present, by the 
cosmic shear and by the compression of the outflowing gas 
into a thin shell, when winds become momentum-driven 
and evolve like a snowplough. In this last case, the den- 
sity of the cooled material accumulating onto the shells, 
which is a mixture of supernova ejecta, shocked interstel- 
lar medium and gas accreted from the ambient medium, 
can become much higher than the density of the ambi- 
ent medium. Shell overdensities 8 are usually in the range 
4/9baryons < 5 < 1000/9baryons • This compression of the gas in 
shells produces and analogous compression of the magnetic 
field lines associated with the outflowing gas, which in turn 
can locally amplify the magnetic fields. We do not attempt 
to model this effect in our simulations. 



4 EXTRAGALACTIC COSMIC RAYS 

Cosmic rays, the other non-thermal component of the in- 
tergalactic and intracluster gas, are relativistic charged par- 
ticles and, as such, they are strongly coupled to magnetic 
fields. Cosmic ray electrons are currently observed through 
the synchrotron radiation that they emit when spiralling 
around the magnetic field lines. Since their cooling time 
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is dominated by the energy losses from Compton scat- 
tering and/or synchrotron emission, energetic cosmic ray 
electrons cool fast, limiting the observability of the pri- 
mary cosmic ray electrons to a short period of their life- 
time. Cosmic ray protons have not yet been unambigu- 
ously detected, but new and up-coming gamma-ray tele- 
scopes, such as the High Energy Gamma Ray Astronomy 
(HEGRA), the Gamma Ray Large Area Space Telescope 
(GLAST), the High Energy Stereoscope System (H.E.S.S.) 
telescope, and the Major Atmospheric Gamma Imaging 
Cherenkov Telescope (MAGIC), are aimed to detect them 
from the gamma rays produced in hadronic interactions 
with the intergalactic and intracluster gas. The relativis- 
tic electrons and positrons generated by hadronic proton- 
proton collisions are one, among others, viable explana- 
tion of the cluster-sized diffuse radio hal oes, observed in 
a significant fraction of galaxy clusters IDcnnison 1980; 
| Vestrandll98llBlasi fc Colafrancescdl 19991: IPolaefc Enfilinl 
2000; IPfrommer fc Enfilinl l200 4a'l . The amount of cosmic 
ray energy required to power those radio haloes is mod- 
erate and comparable to the exp ected magnetic field en- 
ergy (Pfrommcr & Enfilin 2004b). Since cosmic ray pro- 
tons are confined by chaotic magnetic field lines, they are 
likely to accumulate in turbulent regions of the Universe, 
and especially in clusters iVolk^Aliaronian fc^reitschwerdd 
[ 19961: lEnfilin et al]li997l iBerezinskv. Blasi fc Ptuskinlll997t 
Ivolk fc Atovanll2000ri 



4.1 Predictions for Cosmic Rays 

In our own Galaxy, the cosmic r ay energy d ensity is of order 
of the magnetic energy density (|BeckJ|2001J) . In general, the 
cosmic rays that accumulate over cosmic times should be 
relativistic, since sub-relativistic protons can suffer from se- 
vere Coulo mb energy losses in the dense ISM and ICM (e.g. 
lEnfllin et al.lll997ft . If we take this into account, we can ap- 
proximate the equation of state of the cosmic ray gas with a 
polytrope with adiabatic index 7cr = 4/3. This implies that 
the cosmic ray-producing gas suffers exactly the same (rel- 
ative) adiabatic energy losses as the (turbulent) magnetic 
fields. Thus, we can assume that the cosmic ray energy in- 
jected in winds is identical to the magnetic energy density 
and in our wind description we can read any magnetic field 
energy density also as a cosmic ray energy density, that is 
£CR = es ln =B7(8vr). 

Our prescriptions for winds assume that the outflows 
are (initially ) thermally driven, and not cosmic ray-d riven 
JPorfil |2004J . iBreitschwerdt. Voelk fc McKenziel ll99lT) . We 
also implicitly assume that the thermal and the relativis- 
tic components of the plasma outflowing from galaxies are 
not coupled via Alfven waves. If this were to happen, the 
cosmic rays would contribute momentum to the wind and 
an appropriate term Ecn should be included in equation [5] 
This injection of momentum by cosmic rays could in princi- 
ple provide a wind with the necessary kick to escape from a 
galaxy where star formation is not very intense. 

The equipartition of energy between magnetic fields and 
cosmic rays may not be true in all environments. In cosmic 
ray-driven outflows, for example, the cosmic rays leave the 
galaxy without ejecting much gas or magnetic fields. In this 
case, the contribution of winds to the cosmic ray content of 
the ICM and IGM could be much higher than our estimate. 



Hence, our results about the ejection of cosmic rays should 
be treated as lower limits. 

The small energy density in injected cosmic rays can 
also be substantially boosted by the action of structure for- 
mation. It is expected that every gas volume in the IGM 
which ends up in the ICM is passed by shock waves of differ- 
ent strength several times. Many of these shock waves have 
low Mach numbers, implying that they are very inefficient in 
accelerating particles out of the thermal pool. However, any 
pre-existing relativistic particle population is able to gain 
substantia l energy by the passage of even weak shocks (e.g. 
Miniati et al. 2001a b). Furthermore, resonant interactions 
with turbulent plasma waves c an also energise old cosmi c 
ray populations to some degree IIBrunetti et alJl200ll 12004'). 
We do not attempt an estimate of the cosmic ray amplifica- 
tion processes. However, it is possible that also these seed 
cosmic rays may have become an energetic population at 
present days. 



5 CONCLUSIONS 

We have presented new semi-analytic simulations of mag- 
netised galactic winds emerging from star forming galaxies. 
Our set of N-body simulations of structure formation in an 
approximately spherical region of space of diameter 52 h~ l 
Mpc combines a high mass resolution with a region of space 
large enough to allow us to simulate the evolution of the 
winds in their proper cosmological context. We have imple- 
mented new recipes for the transport of magnetic fields into 
the IGM and ICM by outflows on top of the pre-existing 
semi-analytic scheme for the evolution of winds of BSW05, 
which uses the galaxy formation model of ISpringel et alJ 

The assumption that the magnetic fields only weakly 
affect the global evolution of winds is well justified by the 
finding that the total energy of the outflow is several or- 
ders of magnitude higher than its magnetic energy. Here 
we have presented two different scenarios for the transport 
of magnetic fields in winds and for each scenario we have 
analysed two models corresponding to high and low wind 
ef ficiency to blow and trave l far i nto the IGM. As shown 
in iBertone. Stoehr fc Wh ite ( 2005) , the ability of winds to 
travel far into the IGM strongly depends on their mass load- 
ing efficieny: the more mass is accreted from the ambient 
medium, the lower is the probability that winds can blow 
out of galaxies. The mass loading efficiency is parameter- 
ized in the wind model by the constant e. Unfortunately, no 
direct observational evidence is available to constrain this 
paremeter. 

In both scenarios, we have neglected any possible field 
amplification process due to cosmic shear and turbulence in 
clusters. These processes are far too complex to be treated 
correctly in semi-analytic models. In what we called our 
conservative scenario we do not consider any amplification 
mechanism, but we take into account the dilution of the 
magnetic fields due to the adiabatic expansion of the mag- 
netic field lines inside winds. This implies that the results 
of our conservative models should be regarded as lower lim- 
its to any seed magnetisation of the Universe. Conversely, 
in our optimistic scenario we do consider a simple prescrip- 
tion for the amplification of magnetic fields in winds, due to 
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the shear motions of the material inside the wind cavities 
and bubbles. In the two optimistic models, we consider an 
efficient amplification of the magnetic fields by the internal 
shear motions and we set up our prescriptions by assuming 
that the amplification by shear motions counterbalances the 
dilution due to the adiabatic expansion. 

Since cosmic rays can be strongly coupled to turbulent 
magnetic fields, the simulations also allow us to follow the 
seeding of cosmic rays by galactic winds in the IGM and 
ICM. Given our assumption that the magnetic energy den- 
sity equals the cosmic ray energy density, in our simulations 
the cosmic rays are spatially distributed in the same way as 
the magnetic fields. 

We find that the magnetic fields ejected by galaxies with 
stellar masses M* > 10 8 Mq can fill a substantial fraction 
of our simulated volume, producing a mean (seed) magneti- 
sation of order 10 -12 to 10 _s G in the conservative models 
and of order 10 -9 to 10 -7 G in the optimistic models. Mag- 
netic field are not uniformly distributed in space, but rather 
seem to roughly follow the large scale distribution of the 
underlying dark matter density field. 

The ejection of magnetic fields and cosmic rays by galac- 
tic winds is efficient enough to explain the magnetic fields 
and cosmic ray energy densities observed or expected in 
clusters, if additional mechanisms intervene to amplify the 
magnetic fields and to energise the cosmic rays by a further 
factor of 10 to 1000. Amplification mechanisms such as cos- 
mic shear, dynamo effects and others (see Section [3. 31 for a 
discussion) do operate while the fields are ejected by winds 
into the ICM and IGM. It remains to be shown, however, 
whether the cosmic ray protons transported by winds can 
be re-accelerated by intergalactic shocks and turbulence to 
the level expected in galaxy clusters from the existence of 
radio haloes, which might be of hadronic origin. 

Our model does not exclude a priori an additional pri- 
mordial magnetic field. In fact, since at least in the high 
mass loading wind scenario a non-negligible fraction of our 
simulated volume is unaffected by magnetised winds, it may 
in principle be possible to detect primordial magnetic fields 
in those regions which have not been affected by feedback 
effects. 
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tion should both motivate our approach and clearly show 
where its limitations are. 

We first consider the magnetic induction equation: 



— — = V x v x B 
dt 



(Al) 



From this, we easily obtain the evolution equation of the 
magnetic energy density £b = B 2 /(8tt) 



d -B 2 V -v + BB : Vv 

Tt £B ~ 4^ 



(A2) 



where d/dt = d/dt + v • V is the total time (or convective) 
derivative. The magnetic energy stored in the wind is 



Eb = Eb[V w 



d X £B 



Vw (eb)v w 



(A3) 



where (• . -)v„ denotes the average over the wind volume Vw. 
The magnetic field strength of winds provided in the paper 
should be understood in terms of this average. 

We are interested in the time derivative of -Eb[Kv], for 
which we need the variation of the energy density and of 
the volume. Since the latter is the integration volume and 
therefore non-trivially differentiable, one best introduces a 
coordinate system x(t), co- moving with the gas flow, which 
is the solution of the equation x(t) — v(x, t), with the initial 
condition x(to) = xq. This allows a coordinate transforma- 
tion to the static coordinate system a?o, so that the energy 
density in a sub- volume V of Vw evolves accordingly to: 
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(A4) 



The second equation is obtained by transforming back 
to the original coordinate x, while the third equation can 
be derived by using the product rule of determinants and 
explicit calculation. Let us imagine that the wind volume 
can be split into a number iV of sub-volumes V%. within 
each sub-volume the magnetic field can be assumed to be 
statistically homogeneous. We can therefore write 



£b[V w ] = Y. Eb ^ = Yl ^ £B * 



(A5) 



and introduce the rescaled field strength b — B/Bi. The 
time derivative is then 



APPENDIX A: DERIVATION OF THE 
MAGNETIC ENERGY EVOLUTION EQUATION 

In this section we derive the equation for the evolution of 
the magnetic energy in winds. We introduce the set of ap- 
proximations necessary to yield the evolutionary equations 
used in our simulations (Equations 171 and 1151 . This deriva- 



En = 22—EB]yi] + — EB [VKr+i ] 



dt 



dt 



dt 



(A6) 



where the last term accounts for the injection of a new wind 
volume element Vjv+i at the base of the wind. We identify 
this term with Es io = £B in Vin- 

Finally, to derive the simplified Equation|7Jfor the mag- 
netic energy evolution, some further idealised approxima- 
tions are necessary: 
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(i) Vanishing cross correlation between fluid motion 
and magnetic fields. 

(ii) Statistical isotropy of the magnetic fields, which 
implies {bb)vi =|l. 

In order to make the impact of these idealisations more 
transparent, we transform Eq. IA6I so that it can be split 
into an ideal and a non-ideal part, 






Eb — Eb Lii - 



Eb Vw 



(A7) 
where the ideal part (including the left-hand side term and 
the first two right-hand side terms) is identical to Eq. 
The non-ideal part includes three contributions. The first 
term accounts for deviations of the relative expansion of dif- 
ferent regions from the overall expansion rate. This term 
is probably small during the pressure-driven phase of the 
wind expansion, but will lead to a larger inaccuracy of our 
simplified description during the momentum-driven phase, 
when the different regions of the wind expand in a different 
way. The second term accounts for correlations in the ex- 
pansion rate within the sub-volumes with the magnetic field 
strength. For dynamically unimportant fields, we expect this 
term to be small. For dynamically important field strengths, 
we expect the magnetic structures to expand faster than the 
unmagnetized regions, leading to an overestimate of the field 
strength. However, the basic assumption of our description 
is that within the wind volume (which excludes the base of 
the wind) magnetic fields are not dynamically important. 
These two non-ideal contributions vanish if no correlation 
exists between the magnetic fields and the fluid motion both 
on large and small scales. The third non-ideal contribution 
is due to the shearing of the magnetic fields and leads in 
general to a field amplification. This term vanishes if the 
fields are both isotropic and uncorrelated with the fluid mo- 
tion. Otherwise, it should lead to dynamo action, for which 
we account in Eq. 1151 with a very simplified but effective 
description. 
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